Effective level attraction and magnetic flux-induced negative differential conductance 
in two double quantum dot molecules embedded in an Aharonov-Bohm ring 
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We study transport of non-interacting electrons through two quantum dot molecules embedded 
in an Aharonov-Bohm interferometer. The system in equilibrium exhibits bound states in the 
continuum (BIC) and total suppression of transmission. It also shows a magnetic flux-dependent 
effective level attraction and lines of perfect transmission when the intramolecular coupling is weak. 
Out of equilibrium, the current displays two kind of negative differential conductance (NDC) regions, 
which have different origins. One is generated by the usual mechanism of the NDC arising in a double 
quantum dot system. The other is induced by the magnetic flux, and it occurs at small voltages 
and for a well definite range of the intramolecular couplings. We explain this effect in terms of the 
level attraction displayed by the system. 

PACS numbers: 73.21. La; 73.63.Kv; 85.35.Ds; 85.35.Be 
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I. INTRODUCTION 

Advances in experimental techniques at the nanome- 
ter scale have allowed to realize and manipulate quantum 
dots in a controlled waji This has permitted to study in 
these systems a wide spectrum of phenomena 2 , as well as 
to take advantage of some properties of quantum dots to 
build nano-devices, such as rectifiers, amplifiers, lasers, 
and others^^^. For their similitude with atoms, quan- 
tum dots are often viewed as artificial atoms, and two 
or more coupled quantum dots as artificial molecules. 
Multiple quantum dot systems, and in particular the 
double quantum dot (DQD) molecule 7 -, are of particu- 
lar importance, because they are more controllable than 
single quantum dots. Moreover, they admit different 
kind of connections to leads. Initially, most of studies 
considered serial DQDs£, but posteriorly paralle l 9 ' 10 : 11 
and T-shape d 12 ' 13 configurations were also examined, 
emerging quantum interference effects due to the exis- 
tence of different electronic paths. Much attention has 
received the DQD molecule embedded in an Aharonov- 
Bohm ring, where the presence of a magnetic flux adds 
a new tool in the control of the transport properties. 
Aharonov-Bohm (AB) oscillations and Fano effect in a 
DQD molecule embedded in a ring have been observed 
exp er iment all y 9 ' 1 4 : 1 5 . There is much theoretical work 
supporting these experiments and exploring new phe- 
nomena in the DQD embedded in a ring. Fano resonances 
and the magnetic-flux controllability of tr anspor t are ex- 
amined, for example, in Refs. [TtEfillHlillS, both in 
presence and the absence of electronic correlations. The 
existence of bound states in the continuum (BICs) is dis- 
cussed in Refs. [lfj|fl6l The interplay between quantum 
interference and Kondo physics in the parallel-coupled 
DQD has been also explore d 20 : 21 . A closely related sys- 
tem which has received attention is the AB interferome- 
ter with two embedded quantum dot o 22 : 23 . Kubala and 
Konig studied equilibrium transport in this system, find- 
ing an effective flux- dependent level attraction in the lin- 



ear conductance 22 . This effect is caused by the renormal- 
ization of levels by the leads. 

From the applicability point of view, an interesting fea- 
ture exhibited by the transport through quantum dots 
is the negative differential conductance (NDC)2£. NDC 
has been studied in single as well as in double quan- 
tum dot systems, and it has applications in amplifiers 
and oscillators in the microwave, mm- wave and Tera- 
hertz frequency ranges 2 ^. In multilevels quantum dots 
NDC can occur when states have different couplings to 
the leads2§i2L2Si22. In a serial DQD, NDC can be pro- 
duced when the bias breaks the transmission channel ex- 
tended along the system^ .. Other theoretical works on 
generation of NDC in DQD connected in series are Refs. 
I31IJ32IJ33IJ34IJ35I . In a DQD embedded in an Aharonov- 
ring, magnetic-flux-induced NDC was found in the strong 
interdot repulsion regimeA§-. A similar result was found 
by Mourokh and Smirnov in a DQD molecule with three 
terminals 36 -. Recently, it was reported NDC induced by 
the electronic correlation in a side-coupled DQD 37 . 

In this work we study equilibrium and non-equilibrium 
transport through two quantum dot molecules embed- 
ded in an Aharonov-Bohm interferometer. We obtain 
analytical expressions for the transmission and we cal- 
culate numerically the current at zero temperature. In 
equilibrium, the transmission exhibits Fano resonances, 
total reflection, and suppressed peaks as a manifestation 
of bound states in the continuum. Moreover, we find 
for small intramolecular couplings a flux-dependent ef- 
fective level attraction and lines of perfect transmission. 
This effect occurs no matter how weakly coupled arc the 
quantum dots forming each molecule. 

In the non-equilibrium regime, we identify two kind of 
NDC regions in the I-V characteristics, occurring at dif- 
ferent scales and of different origin. One is generated by 
the usual mechanism of the NDC in a double quantum 
dot systems, and it is independent of the magnetic flux. 
The current will increase or decrease with voltage, de- 
pending on whether the voltage makes the levels of the 



2 



different quantum dots become aligned or not aligned. 
The second NDC region is induced by the magnetic ffux, 
and it occurs only for a definite range of intramolecular 
couplings. An abrupt rise of current occurs for small bias 
voltages as consequence of the effective level attraction 
of the hybridized levels produced by the magnetic flux. 
The decrease of current is result of the destruction of this 
effect when the bias voltage is increased. 



II. MODEL 

The system under consideration is shown in Fig. [T] 
Two equal double quantum dot molecules are embed- 
ded in an Aharonov-Bohm ring, which is attached to 
large contacts through one-dimensional leads. Equilib- 
rium transport in a similar configuration with additional 
connections between dots was studied by Li et alM-. The 
left and right contacts are in thermodynamic equilib- 
rium with thermodynamical potentials fiL and fiR, re- 
spectively. The leads are assumed to be ballistic con- 
ductors. We assume that a bias voltage V/e is applied 
between source and drain such that the site energy is V/2 
for the left lead and — V/2 for the right lead. We con- 
sider only one level relevant in each of the quantum dots. 
The system is modeled by a non-interacting Anderson 
Hamiltonian. which can be written as 



H — Hm + H o + Hi, 



(1) 



where Hm describes the dynamics of the isolate 
molecules, 

Hm = Y Y £ A a d ^ 

a=+,- i=A.B 

+ Y ^ d \, a d B, a +tld ] Ba d Aa ), (2) 

<*=+,- 

where £a(b) is the level energy of the left (right) quantum 
dot in the molecule a (a = +, — ); di <a (d\ a ) annihilates 
(creates) an electron in dot i in the molecule a, and t a is 
the intramolecular tunneling hopping. Hq is the Hamil- 
tonian for the noninteracting electrons in the leads, 
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where Cj (c|) is the annihilation (creation) operator of an 
electron in the site i-th of the leads, and v the hopping 
between sites in the leads. The term Hj accounts for the 
tunneling between molecules and leads, 

Hl=- Y ( V aA«C-l + V^^-ldAa) 

a=+,- 

- Y (Vc?d\ a c l+ V**c\d Ba ) 
<*=+,- 
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FIG. 1: Two double quantum dot molecules embedded in 
parallel in an Aharonov-Bohm interferometer 

with Va^ B \ the tunneling coupling connecting the left 
(right) dot of the a-th molecule with the left (right) lead. 
We restrict to the case in which there is not magnetic field 
acting directly on the electrons, so that the situation will 
be identical for the two values of spin. Therefore in what 
follows we omit the spin index. 

In presence of a magnetic flux threading the ring, 
and using gauge invariance, we add the Aharonov-Bohm 
phase (f> — 27r$ / $o around the ring by the replacement 
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(3) 



with $o = h/e the flux quantum. We look for the steady 
states \ipk) of the whole Hamiltonian H. The Hamil- 
tonian describing the leads, Hq, corresponds to a free- 
particle Hamiltonian on a lattice, the eigenfunctions be- 
ing Bloch functions 



(4) 



where \kp) is the momentum eigenstate and \ j) a Wannier 
state localized at the j-th site. The corresponding dis- 
persion relations are e = V/2 — 2v cos k^, for the electrons 
originated in the left contact, and e = —V/2 — 2wcos/cfl 
for those originated in the right contact. The eigenstates 
of the entire Hamiltonian can be written as 
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-E4b> ( 5 ) 

i'=l 



We assume electrons as described by a plane wave inci- 
dent from the far left (right) with unit amplitude, reflec- 
tion amplitude r (r'), and transmission amplitude t (£'). 
Thus, for electrons incident from the left, 



-ik L j 



te 



J < -1 

i>i 



and for electrons incident from the right, 

kR _ [ e - lk M + r ' e ikRj j > 1 

a o ~ U'e- ikz -j, j<-l 



(6) 



(7) 



where the two contributions are independent. Inserting 
Eqs. (O,© and in the Schrodinger equation H\i/)k) = 
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Ek\i>k), we can solve for and b\ a . We are interested 
in the transmission and the current through the system 
for an applied voltage V/e between contacts. We center 
in the the symmetrical configuration, that is, equal left 
and right dot-lead couplings, — ts- Additionally, the 
energies of the quantum dots are £a = Ml and Eb = fiR- 
The current in the leads is given by 



2iev 



where 



C Cj 



= {L,R} k 



(8) 



(9) 



We assume that the voltage drop occurs only between the 
dot A and the dot B, so that the energies of the quantum 
dots are equal to the site energies of the adjacent leads, 
ea = V/2 and £b = —V/2. At zero temperature, the 
states that contribute to the net current are the states 
of the left lead with energies between (j,r = —V/2 and 
/ii = V/2. Evaluating ©-([H]), we arrive to the following 
expression for the current 



2e r v/2 
I(V) = - / T(e) de, 

11 J -V/2 



(10) 



where T(e) is the transmission. We have assumed that 
the Fermi level in equilibrium is equal to 0. 



III. RESULTS 

The transmission probability can be written in the sim- 
ple form 



T(e) 



4t 2 r 2 (Ae) 4 cos 2 0/2 



[(Ae) 4 + r 2 (e - q) 2 } [(Ae) 4 + Y 2 (s + q) 2 } 



(11) 



where (Ae) 2 = e 2 - p 2 , with p 2 = (V/2) 2 + t 2 c and 
q 2 = (V/2) 2 + t 2 c cos 2 0/2), where Y = 4^*^(0) is the 
characteristic line width, with p(0) the density of states 
in the leads at the Fermi level. It is evident from the ex- 
pression that the transmission has a period A0 = 2n 
(A$ = $o)- On the other hand, it is totally suppressed, 
both in and out of equilibrium, when <f> is an odd mul- 
tiple of 7r ($ = n$o/2, n odd). This fully destructive 
interference effect for this value of <f> is expected, since in 
the absence of magnetic flux the upper and lower paths 
are equivalent, then the magnetic flux introduces in the 
wave function a phase —tt/2 along one arm and it/ 2 along 
the other. The total suppression of transmission was 
discussed in an Aharonov-Bohm interferometer with two 
quantum dots, in equilibrium 22 . 



A. Equilibrium transport 

The equilibrium transmission is obtained from Eqs. 
PT|) - (|13p by making V = 0. Let us first analyze the 



action of the magnetic flux in the conductance G 
(2e 2 /h)T(0). This is given by 



G = 



2e 2 4cos 2 (7r$/$ )(i c /r) 2 
~~h [cos 2 (tt$/$ ) + (t c /T)' 



]2 ' 



(12) 



We distinguish in Eq. (fT2")) two different behaviors. When 
t c < Y there is always a value of $ for which G = 2e 2 /h, 
while when t c > Y this never occurs. These behaviors 
are illustrated in Fig. which shows the conductance 
versus the magnetic flux for two different values of the 
t c . For t c = 0.25r, the conductance reaches the max- 
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FIG. 2: (Color online) Conductance versus magnetic flux for 
t c = 0.251" (solid line) and t c = 1.5T (dash line). 

imum G = 2e 2 /h in the interval $ = to $ = $ /2, 
being symmetric around $0/2. Such a maximum occurs 
at e = arccos (i c /r)<i>o/7r. Similar features are found for 
any value of t c < Y, no matter how small it is. This is a 
remarkable result, taking into account that for a molecule 
in series the conductance, proportional a t 2 , gets progres- 
sively smaller as t c decreases^. In contrast, for t c >Y the 
conductance decays monotonously with $ in the same in- 
terval, never reaching the quantum limit. 

The occurrence of maximum conductance for small val- 
ues of t c is associated to an effective level attraction of the 
hybridized levels, similar to the discussed in Ref. |22j in 
an Aharonov-Bohm interferometer with two single-level 
quantum dots. This effect is illustrated in Fig. [31 which 
shows transmission versus energy for t c = 0.25r. In the 
absence of magnetic flux, the transmission behaves qual- 
itatively as that of the DQD in series. In Fig. (dash 
line), t c is small and a single and flat peak is observed. 
The hybridized states are not resolved yet, which is the 
usual for a molecule in series with small t c . However, two 
BICs are occurring, similarly to what happen for a DQD 
embedded in an AB rin g 10 ' 16 . In other words, two of the 
hybridized states are localized and do not participate of 
transmission. This situation repeats whenever $ = n4>o 
(n integer). When the flux is on, the BICs are replaced 
by two Fano resonances, as seen in the rest of Figs. [3] 
The Fano peaks reach T = 1. 

In the sequence from Fig. [3Jl to Fig. [3Ji it can be 
observed the flux-controlled level attraction between hy- 
bridized states. As $ is increased from to $o/2 the 
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FIG. 3: Transmission versus energy for t c = 0.25F and a) 
<E> = (dash line) and $ = 0.1<E> (solid line), b) $ = 0.2$ , 
c) $ = 0.3$o, d) $ = 0.4$ (solid line) and <E> = 0.46$o (dash 
line) . 



Fano peaks get progressively closer to each other until 
overlapping completely, the transmission decaying when 
$ is close to $o/2 (dash line in Fig. [3jl), until vanish- 
ing. The value 2e 2 /h in the conductance shown above is 
reached when the two peaks overlap totally. In contrast, 
for t c — 1.5T the level attraction is not significant, as 
observed in Fig. 2] 
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FIG. 4: Transmission versus energy for t c = 1.5T and a) 
$ = 0, b) $ = 0.03$o, c) $ = 0.1$o, d) $ = 0.2$ (solid 
line) and $ = 0.4$o (dash line). 

It follows from Eq. [11] that perfect transmission takes 
place at energies obeying the following equation 

e 4 - e 2 (2t 2 c - T 2 ) + t 2 c [t 2 c - T 2 cos 2 (0/2)] = 0, (13) 

with two pair of solutions 

ef = ±[A + B 1 ' 2 ] 1 / 2 1 ef = ±[A- B 1 ' 2 ] 1 ' 2 (14) 

with A = t 2 - T 2 /2 and B = Y 2 - 4t 2 c sin 2 (0/2). Fig. 
[5] shows the positions of the T = 1 peaks for different 
values of t c < r/y2- These correspond to ef and e7, 
the only real solutions of Eq. (| 1 3 j) - In the three cases 
the peaks positions shift progressively to the center of 
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FIG. 5: (Color online) Positions of the T = 1 peaks for dif- 
ferent values of t c < T/y2. 



the band as $ increases. The peaks meet at e = at 
<I> = arccos (< c /r)<I>o/7r- The situation slightly changes 
when r/v2 < i c < T, where two new real solutions of 
Eq. (fT3")) arise in the interval arccos (t c /T)&o/ir < $ < 
arccos(l — r 2 /2t 2 )$ /27r, while ef and ej~ become real 
for all $ < arccos (1 - r 2 /2i 2 )<J> /27r. This is illustrated 
in Fig. E for t c = 0.8T and t c = T. The solutions ef 
and ef never meet, but ef and e^ do so at e — 0, being 
responsible of the conductance G — 2e 2 /h for these val- 
ues of t c - When t c > T there are four T = 1 peaks when 
$ < arccos (1 — r 2 /2i 2 )<f> /27r. These peaks never reach 
the center of the band, as shown for t c — 1.5T. 
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FIG. 6: Positions of the maxima of transmission for different 
values of t c > T/\/2. The solid lines correspond to ef and 
and the dash lines to ef and ej. 

We show below that he flux-dependent level attraction 
and the complete overlapping of resonances present for 
small values of t c strongly influence the behavior of the 
non-equilibrium transport. 



B. Non-equilibrium transport 

Let us now consider a voltage V applied between con- 
tacts and let us study the current in the leads. Fig. [7] 
(upper panel) shows the current-voltage characteristics 
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for $ = and different values of intramolecular cou- 
plings. 

In all cases the current-voltage characteristics displays 
a peak, with the corresponding region of NDC. This 
feature occurs analogously to the observed in the serial 
DQD 3 -. The current increases when the bias allows a 
transmission channel exists along the left and right sides 
of the system. If the bias continues to increase, the chan- 
nel is destroyed resulting in the drop of current. On the 
other hand, as occurs for quantum dots connected in sc- 
ries, in general larger values of t c give larger values of the 
current for the same voltage. 

The existence of a magnetic flux produces changes in 
the I-V characteristics which become important when 
t c < r/V2 and $ is within an interval close to $o/2, as 
illustrated in Fig. (lower panel), where $ — 0.46$o- 
We observe sharper current peaks at lower voltages as 
compared to the case $ = 0. The abrupt increase of the 




V/r 

FIG. 7: (Color online) I-V characteristics for $ = (upper 
panel) and $ = 0.46$o (lower panel), and t c = O.ir (solid 
line), t c = 0.25r (dash line), t c = 0.5r (dotted line), t c = 
0.75r (dash-dotted line), t c = F (short dash line). 

current at low voltages for small t c is consequence of the 
level attraction discussed for zero bias. To visualize this 
we have plotted in Fig. [5] the I-V characteristics for a 
fixed t c and <£> (left panel), and the transmission spectra 
associated to the bias voltages indicated in the current 
curve (right panel). As observed, for small bias (cases 
1-3) the transmission keeps large in all the transport re- 
gion (— V/2 < e < V/2), due to the existence of two 
overlapped resonances close to each other. Larger bias 
voltages make the heights of the resonances fall, so that 
the transmission in all the window of transport becomes 
smaller, occurring the observed decrease in the current. 
A further insight of this is obtained through the density 
of states of the left and right quantum dots. Fig. [5]shows 
the left and right quantum dots DOS for the same param- 
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FIG. 8: (Color online) I-V characteristics for t c = O.ir and 
$ = 0.46$o (left panel). Transmission spectrum for different 
bias voltages for the same parameters (right panel). 

eters of Fig. [8l In equilibrium (V = 0) a molecular state 
is formed. For the cases 2 and 3 the coherence is still 
preserved but for higher voltages (cases 4-5) the physical 
picture changes. In these cases, the coherence between 
dots is lost, the electron is localized at the left quantum 
dots and the molecular bridge is broken. 
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FIG. 9: (Color online) Densities of states of the electrons 
coming from the left at the quantum dots A (left panel) and 
B (right panel), for t c = O.ir and $ = 0.46$ . 

Fig. [TUlshows the current versus bias voltage and mag- 
netic flux for t c = O.ir. It can be identified clearly the 
two regions of maximal current at different scales. In the 
upper panel, there is a broad peak centered in $ = 
and V ~ I.2r. The behavior of the current versus volt- 
age in this region was already discussed for zero flux in 
Fig. [7] The current maximum in this case decreases 
monotonously with the magnetic flux, remaining its po- 
sition almost unchanged. In the same panel it is high- 
lighted a region close to $o/2 and small voltages, which 
is plotted in the lower panel. The observed peak corre- 
sponds to the enhancement of current with the magnetic 
flux taking place for weak intramolecular couplings, de- 
scribed in Fig. lower panel. The maximum of current 
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now occurs around $ = 0.44$ and V — 0.0291". Also, 
this peak is clearly defined when t c < T/y/2, when the 
flux produces a total level attraction in equilibrium and 
it is not present when t c ~ T or larger. It is important 




to note that the latter feature does not exist if the mag- 
netic flux is absent, so that in this case we can properly 
talk of magnetic flux-induced NDC. Similar results are 
discussed in a parallel DQD molecule embedded in an 
Aharonov-Bohm ringi^ and in a molecule in a three ter- 
minals configuration^, in both cases NDC occured in the 
strong interdot repulsion limit. 

We expect that the above picture remains valid even 
if the electron-electron interaction is taken into account. 
In fact, in embedded QD arrays, the main effect of the 
electron-electron interaction is to shift and to split the 
resonance positions 7 7 . This occurs because the on-site 
Coulomb repulsion energy U introduces a renormaliza- 
tion of the site energies. In analogy with QD arrays in 
series, we expect that depending on the relation between 
the interdot coupling and the on-site Coulomb interaction 
different regimes arise. For t c /U -C 1, the resonances 
and antircsonances would split into two distinct mini- 
bands separated by the on-site Coulomb energy, while 
for t c /U 1, the resonances and antiresonances would 
occur in pairs. We think that the above behavior would 



not break the negative differential conductance. A work 
in this direction is under progress. 



IV. SUMMARY 

We studied the transmission and the I-V characteris- 
tics for two double quantum dot molecules embedded in 
an Aharonov-Bohm ring. We showed that for t c < T, 
the magnetic flux can be used to control totally the con- 
ductance, allowing this to take any value between and 
2e 2 /h. When t c < T/\/2 the flux produces an effective 
level attraction and lines of perfect transmission, allowing 
the levels to meet at the center of the band at a deter- 
minate value of the flux. On the other hand, the system 
displays BICs when $ = n$o ( n integer), and total sup- 
pression of transmission when $ = n<f>o/2 n odd. In the 
non-equilibrium regime, we identify two kind of NDC re- 
gions in the I-V characteristics, occurring at different 
scales and of different origin. A first current peak exists 
at voltages of the order of the characteristic linewidth T, 
and it is independent of the magnetic flux. The role of 
the flux in this case is to control the height of the overall 
current. The drop of current with the increase of bias 
has analogous explanation to the NDC region in a se- 
rial DQD. A second peak in the I-V characteristics takes 
place only when t c < T, at voltages V ~ T/10 and it 
is strongly dependent on the magnetic flux. In fact, it 
does not exist if the flux is absent. The current suffers 
an abrupt rise for small bias voltages, as consequence of 
an effective level attraction of the hybridized levels pro- 
duced by the flux. The decrease of current is result of 
the destruction of this effect when the bias is increased. 
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